Acoustic Faraday effect in Tb 3 Ga 5 Oi2 
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The transverse acoustic wave propagating along the [100] axis of the cubic Tb3GasOi2 (acoustic 
C44 mode) is doubly degenerate. A magnetic field applied in the direction of propagation lifts this 
degeneracy and leads to the rotation of the polarization vector - the magneto-acoustic Faraday 
rotation. Here, we report on the observation and analysis of the magneto-acoustic Faraday-effect 
in Tb3GasOi2 in static and pulsed magnetic fields. We present also a theoretical model based 
on magnetoelastic coupling of 4/ electrons to both, acoustic and optical phonons and an effective 
coupling between them. This model explains the observed linear frequency dependence of the 
Faraday rotation angle. 



PACS numbers: 62.65. +k, 72.55. +s, 73.50.Rb 

I. INTRODUCTION 

TGG is a dielectric material with a cubic garnet struc- 
ture. The orthorhombic local symmetry for Tb 3+ (J = 6) 
leads to pronounced crystal electric field (CEF) effects 
and a Curie-type magnetic susceptibility. At Tn = 
0.35 K an antiferromagnetic transition appears^— 

Remarkable symmetry experiment, the so-called 
phonon Hall effect, was recently performed and anal- 
ysed in Tb3GasOi2 (TGG)^- 7 - : in analogy to the classical 
Hall effect in conducting materials, the appearance of a 
thermal gradient in the direction perpendicular to both, 
the applied magnetic field and the thermal flux, was ob- 
served. This observation suggests existence of magneto- 
acoustic phenomena in TGG. One of them, the acoustic 
Farady effect, has been observed and is reported here (see 
also Ref. §). 

The intriguing magneto-elastic properties of TGG can 
be described using the magneto-elastic interaction of the 
form 

i? = X>e r Or, (1) 
r 

which represents the basic interaction mechanism of the 
phonon modes with the orbital moments of the Tb 3+ 
ions, where T denotes the symmetry label, er - the strain 
components, Or - the quadrupolar operators of the Tb 
ion, and gr - the magneto-elastic coupling constant. 

The elastic constants en, (en— ci2)/2, and C44 in TGG 
exhibit strong anomalies below 100 K<£ These anomalies 
are clear evidence for the CEF magneto-elastic interac- 
tion described by Eq. (1). The large coupling constant 
for the C44 mode g(T§) points to the possibility for inves- 
tigating special magneto-acoustic effects, among others, 
the acoustic Faraday effect. This effect was also observed 
previously in ferrimagnetic YsFesO^ (YIG), in antifer- 
romagnets, and paramagnets. For a review see Ref. [Ioj |. 



In analogy to the optical Faraday effect, for the 
magneto-acoustic Faraday effect a transverse acoustic 
wave propagating along a fourfold cubic axis has two, left 
and right circularly polarized, components. In an applied 
magnetic field these components have different velocities. 
This leads to a rotation of the linear polarization vector 
as a function of magnetic field, seen as oscillations in 
amplitude of the acoustic signal. For the optical Fara- 
day effect, the rotation angle per unit length is given by 
$>/L = VB, with V being the Verdet constant. We will 
show that for the magneto- acoustic Faraday effect, $/L 
can be more complicated. 

II. EXPERIMENTAL RESULTS 

The experiments have been performed on a TGG sin- 
gle crystal oriented for propagating the sound wave with 
wave vector k along the [100]-direction. The sample 
length along the direction of the sound-wave propagation 
was 4.005 mm. The ultrasonic attenuation and velocity 
have been measured with a set-up as described at great 
length in Ref. [lj| [H|. LiNbOa transducers have been 
used with polarization u along the [010] direction. Fields 
up to 20 T have been provided by a commercial supercon- 
ducting magnet system. Experiments beyond that field 
have been performed in pulsed magnets. Both systems 
have been equipped with 4 He-flow cryostats. 

Static fields 

In static magnetic fields, we have performed measure- 
ments for frequencies 97, 173, 177, 215, and 337 MHz 
at constant temperature of 1.4 K. Figure [1] shows typi- 
cal patterns of the magneto-acoustic Faraday oscillations 
observed for TGG at different frequencies. The ampli- 
tude is plotted versus field. For the fields above about 
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FIG. 1: (Color online) Amplitude oscillations of the acoustic 
C44 mode versus static magnetic field in Faraday geometry 
B || k || [100], u || [010] measured at frequency of 97, 173 
and 215 MHz and constant temperature of 1.4 K. The arrows 
indicate the field sweep direction. 



13 T the transducer detects subsequent maxima for po- 
larization angles $ = nn and subsequent minima for 
<3> = (2n+ l)7r/2 in intensity or amplitude. Near 20 T the 
oscillations cease. In pulsed-ficld experiments we have 
observed that the oscillations reappear above 20 T (see 
next subsection). 

In Fig. [2] we plot the rotation angle per unit length 
and frequency $/L/o, where L is the acoustic path length 
for a particular echo (sample length in our case) and /o 
is the measurements frequency. Presuming a linear fre- 
quency dependence, the $/L data for different frequen- 
cies have been divided by the measurements frequency Jo- 
lt is seen that indeed the linear frequency dependence is 
fulfilled to a high degree of accuracy. Such a linear de- 
pendence was observed before in the paramagnetic state 
of CeAl2.— This result is very astonishing because ear- 
lier theoretical treatments of the magneto-acoustic Fara- 
day effect, for which only the magneto-elastic coupling 
of Tb 4/ electrons with acoustic phonons has been taken 
into account, show always a quadratic frequency depen- 
dence for B 3> 27r/o/7 with 7 being the gyromagnctic 
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FIG. 2: (Color online) Rotation angle per unit length and 
frequency &/Lfo as a function of magnetic field for frequen- 
cies 97, 173, 177, 215 and 337 I 
amplitude maxima from Fig. Q] 
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This condition is easily obeyed for B > 10 T 



FIG. 3: (Color online) The relative velocity change for the 
acoustic C44 mode as a function of magnetic field for T — 1.4 
K, B || k || [100], polarization vector u || [010]. Symbols - 
experimental data, dashed line - theoretical results. 



and frequencies in the 100 MHz range. For 100 MHz, the 
cyclic frequency uj = 2nfo, and the ratio w/7 = 3.7 mT. 
This discrepancy was resolved recently^ and is discussed 
later in the theory section. 

In addition, we measured the field dependence of the 
velocity of the C44 mode at fixed temperature of 1.4 
K (Fig. [3]). This mode exhibits a strong softening 
with a pronounced minimum at about 18 T, i.e., at the 
field where the lowest level of the quasi-triplet mode 
crosses the upper doublet mode (see Section 3). Here 
the maximum softening of the sound velocity amounts to 
Av/v = 6%. 

We describe Av/vq by calculating the strain susccpti- 
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bility in the magnetic field using the same quadrupolar 
operator Or = O xy + O.6O2 as for the temperature de- 
pendence of C44.— »i£ We use the formula 



cr = c r (B = 0) - -g^Nxv 



with measured elastic constant c r 



P v o 



2 - 9.67 x 



10 10 J/m 3 . Here p is the mass density. Further, N = 
1.28 x 10 22 cm -3 is the number of Tb-ions per cm 3 , and 
Xv is the strain susceptibility, which is defined in analogy 
to the magnetic susceptibility as xr = d(Or)/der- 

There is a rough agreement of the overall shape of 
the calculated Av/vq (dashed line) with the experimen- 
tal curve (solid line), though the calculated minimum 
at 21 T is at slightly higher field than the experimen- 
tally observed minimum in static fields. As for the cou- 
pling constant we surprisingly find a value of gr ~ 21 K, 
which is roughly a factor of two smaller than the coupling 
gr = 45.5 K determined from a fit of C44(T).— 

We tentatively attribute this reduction to the fact, 
that the CEF levels display a nonzero (possibly field- 
dependent) linewidth, which is not included in our cal- 
culation. Also we expect that a careful averaging over 
all incquivalent ions, which will force us to incorporate 
new strain coupling constants, will modify the picture to 
a certain degree. In the above calculation we took this 
averaging into account by an overall factor |, because 
we found Xr(T) to be strongly suppressed for one third 
of the ions. Details of the calculation and possible ef- 
fects of the averaging procedure will be presented in a 
forthcoming publication^ 



Pulsed fields 

To confirm the reappearance of the acoustic Faraday 
oscillations above 20 T we further performed pulsed-field 
experiments. In Figs. U and [51 we show the field de- 
pendence of the sound amplitude oscillations and sound 
velocity taken at 93 MHz. The data demonstrate some 
hysteresis. The acoustic Faraday oscillations do occur 
both below and above 20 T. 

However, some puzzling phenomena appear in pulsed 
fields. The oscillations in Fig. @] can be observed in 
the complete field range with the exception of a damped 
region around 20 T. This is in sharp contrast to the 
static-field measurements shown in Fig. [TJ where the 
oscillation period is very large for B < 13 T. In addition, 
the decrease of the corresponding sound velocity is 30 % 
stronger in pulsed fields (Fig. [5j than in static fields 
(Fig. [3]). These differences between static- and pulsed- 
field data possibly are due to a nonequilibrium situation 
during the pulsed-field experiment. The time dependence 
of B is given in the inset of Fig. 2J For the static field 
measurements we swept the field at a rate of 75 mT/min. 

Under quasi-adiabatic conditions, a huge magne- 
tocaloric effect (MCE) can occur in pulsed fields. In 




FIG. 4: (Color online) Amplitude oscillations of the acoustic 
C44 mode versus pulsed magnetic field in Faraday geometry 
B || k || [100], u || [010] measured at frequency of 93 MHz and 
temperature of 1.4 K. The arrows indicate the field sweep 
direction. Inset: time profile of pulsed magnetic field. 
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FIG. 5: (Color online) The relative velocity change for the 
acoustic C44 mode as a function of magnetic field for T — 1.4 
K, B || k || [100], polarization vector u || [010]. Arrows 
indicate direction of the field sweep. 



TGG, the MCE effect was investigated theoretically for 
fields applied in different crystallographic directions and 
experimentally for B || [110] j 17 i 18 We performed calcula- 
tions of the MCE as in Ref. [lU , averaging over the CEF 
levels of the three Tb 3 + ions in incquivalent positions. In 
Fig. [6] we present our results for different starting tem- 
peratures and for magnetic field applied along the [100] 
direction. From our calculations we expect nonmono- 
tonic temperature changes with magnetic field, as seen 
in Fig. HI The MCE, together with thermal relaxation 
to the bath may lead to complicated hysteretical heat- 
ing and subsequent cooling effects. This could at least 
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FIG. 6: (Color online) Temperature change due to the MCE 
for B || [100] calculated for different initial temperatures of 
To = 0.3, 1, 2, 4 and 10 K. 
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FIG. 7: (Color online) Rotation angle $, extracted from data 
in Figs. [4] and [TJ for pulsed and static fields, respectively. Ul- 
trasound frequency is 93 MHz in pulsed-field measurements; 
95 MHz - in static-field measurements. 



partially explain the increased amplitudes of Faraday os- 
cillations in the descending held. A quantitative analysis 
of these effects is challenging and stays outside the scope 
of this work. 

From the acoustic amplitude oscillations of the pulsed 
field data the rotational phase versus magnetic field can 
be extracted. This has been done in Fig. [7] for oscilla- 
tions in the field range B < 20 T and B > 20 T. This fig- 
ure qualitatively agrees with the theoretical predictions, 
shown in Fig. [8] 

The obtained results for the magneto-acoustic Faraday 
oscillations and the sound velocity-change are a good il- 
lustration that various side effects have to be taken into 
account in pulsed-field experiments. 



FIG. 8: The calculated normalised Faraday rotation angle 
$/$o('I > o = ^k(m,QEg) 2 uj ; k = to/v) as function of mag- 
netic field for T < A t , B || k || [001] ± u. Full line: using 
bare a-o coupling E; broken line: using renormalized Et with 
2mQ§ 2 /E = 0.05. niQ is the quadrupolar (O xz ,O yz ) singlet- 
triplet matrix element for A* excitations. Finite linewidths 
for uj n and Af have been used. 



III. THEORY 

In TGG, the cubic crystal leads to fourfold symmetry 
axes [001] etc. Therefore, doubly degenerate transverse 
acoustic (a) and optical (o) lattice vibration (E) modes 
exist for a wave vector k parallel to these axes. Equiv- 
alcntly, they may be presented by left (L) and right (R) 
circularly polarized modes, which arc complex conjugates 
and degenerate due to time-reversal symmetry. Applica- 
tion of a magnetic field breaks this symmetry and leads 
to different wave numbers for L and R modes of a given 
frequency uj, provided that there are magnetic degrees of 
freedom, which couple sufficiently strong to the lattice vi- 
brations. The former may be collective spin waves as well 
as localized CEF excitations. In the case of spin waves in 
a ferro (fcrri-) magnet the simplest theoretical treatment 
of the magneto-acoustic Faraday effect takes into account 
the coupled equation of motion for the classical magnetic 
moments M and the acoustic displacements u. Then one 
obtains a Faraday-rotation angle per unit length given by 
the difference of the L and R wave numbers a a 10 ' 19 



$ 1 1 

— = — (kr — kfi) = —sk , „ 
L 2 V ' 2 \uj 2 



(2) 



Here, LO m = jB e g and s = jb 2 / (v 2 Mq) , where i? c ff is the 
effective molecular field, Mo stands for the magnetic mo- 
ment in the propagation direction and 7, b and v are the 
gyromagnetic ratio, the magnetoelastic coupling constant 
and the sound velocity, respectively. For large enough 



fields u> m 3> LJ and then <1>/L oc uj , since k 



Calcu- 



lations for S = 1/2 paramagnetsi^ and multilevel param- 
agnetic CEF systems^ also give a dependence $/Lau 2 . 
Thus, existing theories so far consistently predicted a dif- 
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FIG. 9: (Color online) Lowest CEF energy levels for field 
applied in [100] direction. Notice the crossing of the lowest 
triplet level with the upper doublet level. 



ferent frequency dependence than the linear frequency 
dependence of $/L observed either in (paramagnetic) 
CeAl?^ or now in TGG. The Faraday rotation in mag- 
netically ordered compounds such as YI C 20 i 21 and para- 
magnets with magnetic impurities^ 2 - was only measured 
for a single frequency. We suspect that also in these cases 
the discrepancy to the quadratic frequency dependence 
of Eq. exists. In this equation, the prefactor k = ^f- is 
unavoidable because it counts the number of sound-wave 
periods per unit length. The problem lies with the factor 
uj in the numerator of Eq. ([5]) which stems from the non- 
diagonal (xy) and purely imaginary component of the 
dynamic magnetic susceptibility that has to vanish for 
uj — > 0, i.e., the splitting of the acoustic L,R modes has to 
approach zero for k — > 0. Thus, it seems that any mecha- 
nism involving only acoustic phonons and magnetic exci- 
tations will lead to the $/L oc u 2 frequency dependence 
in contrast to experiment. A solution for this discrepancy 
was suggested recently^ In non-Bravais lattices such 
as TGG (and also CCAI2) k -> optical phonons may 
have a direct influence on sound waves due to a coupling 
of elastic long-wave length strains and internal (optical) 
displacements, provided they have the same symmetry. 
For TGG, there are many optical phonons of the right 
(doubly degenerate) E symmetry 2 ^, that can couple to 
the E-type transverse acoustic strain components^ The 
acoustic phonons have a direct magnetoelastic coupling 
to CEF transitions (Eq. ([1])), which again leads to a Fara- 
day rotation with quadratic frequency dependence of the 
generalized form of Eq. Furthermore, the optical 

phonons also couple to the CEF excitations. This leads 
to a splitting of optical E-type (transverse) phonons with 
bare frequency ui into circular polarized modes with split 
frequencies uj^ 2 . Via the effective coupling of acoustic 
and optical modes this creates an indirect contribution 



L ~ 2 kEt 2 



~-2 
Wo 



A.) 

r,+ 2 J ' 

Ulo ' 



(3) 



with circular optical phonon frequencies and effective a-o 
coupling E t given by 



^ 2 



^(i+f ((o xz o xz )y Wo ±~9 2 ((o xz o yz ))Z o ) 

E t = E + g((d xz d xz ))' (4) 



Here, g a ,g ,g = (ffaSo) 2 are acoustic, optical, and 
mixed magnetoelastic coupling constants and E is the 
bare a-o coupling of phonons. Furthermore, the dou- 
ble brackets denote the dynamical local susceptibilities 
(prime: real part, double prime: imaginary part, sub- 
script: frequency) of 4f quadrupolar-moment operators 
O xz = (J x J.z + JzJx) and O yz = {J y J z + J z J y ) which 
couple to x, y polarized acoustic and optical phonons. 

For k — > in Eq. ([3]), the splitting of optical phonon 
modes stays finite and, therefore, (k = ui/v) we have 
a linear frequency dependence <I>/L ~ uj for the indirect 
Faraday-rotation angle induced via a-o phonon coupling. 
It has been arguecU^ that for reasonable parameters this 
term dominates the quadratic term in Eq. 

For the explicit calculation of the contribution linear 
in lo we need the above quadrupolar susceptibilities for 
the CEF level scheme of TGG in a magnetic field. Due 
to the low orthorhombic symmetry the J = 6 manifold 
splits into 13 singlets 2 ^ at three incquivalent Tb 3+ sites. 
The Zccman splitting of the lowest levels of the most 
relevant inequivalent Tb 3+ ion for B || [100] is shown 
in Fig. [9l Apparently, some quasi-degeneracies remain 
for zero field, e.g. an excited magnetic quasi-triplet at 
At ~ 57 K. Therefore, the low-lying levels may be ap- 
proximately described by a simplified CEF model with 
cubic symmetry as proposed in Ref. 0. It is used here 
for the calculation of the quadrupolar susceptibilities in 
Eq. (|4|). Their explicit expressions for the cubic doublet 
(ground state) - triplet (excited state) model are given 
in Ref. [l5j]. The result of a model calculation for the 
Faraday rotation using Eq. (TJJ and Eq. ^ is shown in 
Fig. [8] 

It is qualitatively very similar to the result of the 
pulscd-field experiment in Fig. [7] The main CEF tran- 
sitions contributing to the Faraday rotation in the cu- 
bic model are those between the ground state and the 
two linearly split (quasi-triplet) excitations with ener- 
gies A^ indicated in Fig. |H1 While A^" contributes to 
the effective a-o interaction E t , A^T leads mainly to the 
splitting of optical uj^ modes. The Faraday rotation 
develops a resonance-like behavior due to two effects: 
The ground state - triplet level crossing (A^(B) = 0) 
and the resonance with the optical phonon frequency 
(A t "(_B) = lo ) which happen at approximately the same 
field of 17 - 20 T. 
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IV. SUMMARY 

In summary, we have observed and analyzed the 
magneto-acoustic Faraday effect in TGG in magnetic 
fields up to 42 T. A strong damping of the Faraday os- 
cillation amplitude has been found in the region of the 
CEF level crossing at about 20 T. We have also shown 
that the rotation angle of the shear-wave polarization 
vector is strictly linear in the sound-wave frequency. The 
strong softening of the acoustic C44 mode as a function of 
magnetic field has been detected and reproduced quan- 
titatively in a calculation taking into account magneto- 
clastic interaction. For the acoustic Faraday rotation our 
theoretical model is based on magnctoclastic coupling of 



4/ electrons to both, acoustic and optical phonons and an 
effective interaction between them. The finite splitting of 
optical phonons for long wavelengths and the acoustic-to- 
optical phonon coupling results in a contribution to the 
Faraday-rotation angle which is linear in the sound-wave 
frequency^ as it was observed in the experiments. 
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